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We revisit the protocols to create maximally entangled states between two Josephson junction 
(JJ) charge phase qubits coupled to a microwave field in a cavity as a quantum data bus. We devote 
to analyze a novel mechanism of quantum decoherence due to the adiabatic entanglement between 
qubits and the data bus, the off-resonance microwave field. We show that even through the variable 
of the data bus can be adiabatically eliminated, the entanglement between the qubits and data 
bus remains and can decoher the superposition of two-particle state. Fortunately we can construct 
a decoherence-free subspace of two-dimension to against this adiabatic decoherence. To carry out 
the analytic study for this decoherence problem, we develop Frohlich transformation to re-derive 
the effective Hamiltonian of these system, which is equivalent to that obtained from the adiabatic 
elimination approach . 

PACS numbers: 73.21. La,03.65.-w, 03.67.Ca, 76.70.Cr 



I. INTRODUCTION 

As a useful quantum resource, entanglement can not 
only be used to test fundamental principles in quantum 
mechanics, such as Bell's inequalities, but also play a 
central role in quantum information processing including 
quantum computation, quantum teleportation and quan- 
tum cryptography. Therefore how to create a stable and 
controllable entangled state in a quantum bits (qubit) 
system is very important for quantum information pro- 
tocols 

A number of protocols have been proposed to produce 
quantum entanglement in different qubit systems, such 
as NMR, polarization photon, quantum dots, Josephson 
junction. Due to the prompt progresses in preparing 
various solid state qubits, these schemes become very 
promising to realize the practical quantum computing. 
Actually according to the DiVincenzo criteria the cou- 
plings JJ qubits 2] for quantum computation, the solid 
system is one of the best candidates for quantum compu- 
tation, since qubit should be scalable, controllable and 
with longer decoherence time. Actually it seems diffi- 
cult to fulfill all the requirements by quantum informa- 
tion processing. Recently several groups have demon- 
strated the macroscopic quantum coherence of Joseph- 
son junction (JJ)qubits with long decoherence time in 
experiments 0] • 

Quantum entanglement plays the central role in inte- 
grating multi-qubit to form a scalable quantum comput- 
ing. We notice that, in most of the protocols to produce 
such J J qubit entanglement, and correspondingly to carry 



'Electronic address: lichong@itp.ac.cn 



out two qubit logic gate operations, each qubit interacts 
with a common quantum object as a data bus, which may 
be an electromagnetic cavity field, a quantum transmis- 
sion line coplanar cavity or an nano-mechanical resonator 
[g| [j| . If the characterized frequency of the quantum data 
bus is off- resonate to the energy spacing of the qubit, the 
degree of freedom of the quantum data bus and the vari- 
ables of the quantum object can be separated adiabati- 
cally form that of two qubit system. Then the induced 
inter-qubit interactions can create an efficient quantum 
entanglement of two qubits. 

However, as we have investigated [Tfij . there usually ex- 
ists quantum entanglement between the states of data 
bus and those of the two qubit system even after re- 
moving the data bus. This adiabatic quantum entan- 
glement has been studied according to the generalized 
Born-Oppenheimer (BO) approximation^^ where the 
slow variables can be driven by different effective poten- 
tials provided by the fast internal states and then the 
entanglement between fast and slow variables forms. Re- 
cently, Averin et al similarly considered the adiabatic 
entanglement of two JJ charge qubits. In this investiga- 
tion, two JJ charge qubits arc assumed to be coupled with 
a large junction which works as a faster data bus. With 
the BO adiabatic approximation, the energy of the low- 
est band of the latter junction can be considered as the 
effective interaction between the two JJ charge qubits. 
But the quantum decoherence induced by the adiabatic 
entanglement has not been considered here though it may 
occur in the case with higher excitation of large junction. 

In this paper, we are also specific to the J J qubit sys- 
tem coupling the cavity and show that the adiabatic en- 
tanglement may cause the extra errors of the logic gate 
operation for this two qubit system with high-excitation. 
We will consider decoherence of the JJ qubit caused by 
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FIG. 1: SQUID SI and SQUID S2 in a cavity coupled to a 
microwave field. 



the the thermal excitation of large junction through this 
adiabatic entanglement mechanisim. Actually, without 
considering thermal excitation, we are not clear if the 
created entanglement between two JJ qubits are stable 
since it can be produced according to an effective Hamil- 
tonian, which is obtained in usual by "ignoring" inter- 
meddle variables of data bus ^1 • 

To carry out a totally analytic study, we utilize the gen- 
eralized Frohlich transformation to re- derive the effective 
Hamiltonian of this system. In this way we can study in 
details this novel decoherence phenomenon for the en- 
tanglement of two JJ-qubits. There exist four entangled 
states for two JJ qubit system, including two maximally 
entangled states that can be obtained by controllable the 
micro-wave field. Though the superposition of some two 
qubit states can decoher due to the adiabatic entangle- 
ment, there exist a decoherence-free subspace, against to 
the decoherence induced by the adiabatic separation pro- 
cess. Therefore, we found that only two of four maximal 
entangled states are stable in this scheme. 



II. THE MODEL OF TWO JJ QUBITS IN 
CAVITY 

Without loss of generality, we investigate a simplified 
model, which consisting of two J J qubits in a cavity with 
a single mode micro- wave fields (FIG. I). The Hamilto- 
nian of the coupled system H can be described as a sum 
of that of the junctions, the cavity field and a interaction 
term between the cavity and the junction @|9j, i.e, 

H = tiuia^a + 4Eci (ni — n g i) 2 — Eji (&) cosifi 

+4E C2 (n 2 - n g2 ) 2 - E J2 ($) cos ip 2 , (I) 

where 

E c = e 2 /2(C g + 2Cj) (2) 

is the single-particle charging energy of the island, Cj the 
capacitance of the junction, C g the capacitance of gate 
and ip i the phase difference between points on the oppo- 
site sides of the i-th junction. The Josephson coupling 
energy 



Ej ($) = 2E J0 cos 



2tt$ 



(3) 



depends on the the total flux $ and the maximal coupling 
energy Ej = I c ^. Here, I c is the critical current of the 
junction, and $o the total flux and flux quanta. 

When a nonclassical microwave field with the vector 
potential 



A (r) = ~u a (r) a + ~u* x (r) a* ■ 



(4) 



is applied, where and a are the creation and annihi- 
lation operators of the cavity fields, the total flux $ is 
divided into two part 

<F = $ e + $/. (5) 

$6 is static magnetic flux through the SQUIDs and 

$ / = |$ A |(e- l9 a + e J V), (6) 

the microwave-filed-induced flux through the SQUIDs 
where 



$A = 



~u a (r) • d I . 



(7) 



We take E C i = E C 2 and Ej\ (<&) = E J2 ($) and then 
the Hamiltonian of the system becomes 



( 7T<f> e + <&f \ , i 
- IE jo COS ( ^— ^ ] (a! 



(8) 



where the quasi-spin operators a x , a y , and a z are defined 
with respect to the the states |0) and |1) of no (one) 
excess cooper pair on the island. To form a qubit or 
a two-level system, one need to tune the gate voltage 
V g so that n g is approximately a half-integer. In this 
case the charge eigen-states are |0) and |1). We assume 
I^aI <S $0) and focus on the charging regime Ec ^S> Ej. 
Then, the Hamiltonian can be approximated as 



H = H + H! 



Ho = tiwj {a\ + fi^) + hoja'a 
Hi = g(a + a t ) (cr* + a'l) 



(9) 



where 



hujj = 2E C 



3t^-(2n + l) 



and the coupling constant between qubit and the cavity 
field is 



2tt V 2/i J s 4>o 
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In practice we take the volume of the cavity and 
the wavelength of microwave respectively as ~lcm 3 and 
~lcm, the dimension of the Josephson junction as ~l/zm, 
the critical current of the junction as / C ~10~ 5 A. Due 

to Eq(4), we have — <§; 1 which means Hi <C 

Ho. So we can perform perturbation theory repre- 
sented by a generalized Frohlich transformation[lj| on 
the Hamiltonian(2). Then we can obtain the effective 
Hamiltonian of two JJ qubit by removing the variables 
of the microwave field approximately. 

III. THE EFFECTIVE HAMILTONIAN FROM 
THE GENERALIZED FROLICH 
TRANSFORMATION 

In its original approach for superconductivity BCS the- 
ory, the Frohlich transformation 0] is utilized to get 
the effective Hamiltonian for electron-electron interac- 
tion from electron-phonon interaction. In general we can 
consider a interaction system described by a sum of free 
Hamiltonian and interaction Hamiltonian, 

H = H Q + Hi. (11) 

Comparing with the free part Ho,the interaction part Hi 
can be regard as a perturbation. Let us define an anti- 
Hermitian operator S, and a corresponding unitary op- 
erator C/t = cxp{— S}. We perform an unitary transfor- 
mation on the Hamiltonian(ll) by this unitary operator, 
and then get the equivalent Hamiltonian as 

H = U ] HU 

= H o+Y,£rm [s ' [ ''' [s ' [s ' Hl]] ''' ]] - (12) 

n— 1 ^ ' * 

n 

Since the unitary transformation U is time - 
independent, Hamiltonians (11) and (12) describe the 
same physical process. We can take both the interaction 
Hi and operator S in the first order terms in the right 
hand side. At the same time, we require the operator S 
to satisfy the following condition 



where \m) and E m are the eigenvectors and eigenval- 
ues of Hq respectively. The transformation, by which 
one can draw out effective Hamiltonian (14) from the 
Hamiltonian(ll), is the so-called general Frohlich trans- 
formation. It has been proved in Rcf 13] that this gen- 
eralized Frohlich transformation is just equivalent to the 
second-order perturbative theory. 

Now we use the above approximation method to de- 
rive the effective Hamiltonian for the two-JJ qubit en- 
tanglement. Under the condition g <C huj, we explicitly 
construct the anti-Hermitian operator S of Eq. (3) as fol- 
lowing 

S = £{A + (a-ai){al + al) 

+iA_ (a + af) (4 + 4)} ( 16 ) 
where the coefficients 

± i^hu) — 2Hu)j hu + 2ftwj ) ^ ^ 

Using the above explicit expression for anti-Hermitian 
operator S, we can finish the generalized Frohlich trans- 
formation and then obtain effective Hamiltonian obvi- 
ously. 

HcS = ^ (al + a 2 ) - (a + a^f {a\ + a*) 

2 2 

+Wa-|-A + -|-AX^- (18) 

If the micro-wave field is very weak, we can discard the 
second terms of a 2 and a* 2 in the effective Hamiltonian 
under the rotating wave approximation. Then the effec- 
tive Hamiltonian Eq.(18) reads 

H cff = (fiwj + g 2 A-a)a) {a\ + a 2 ) 

2 

+hu)a?a + y A + ^cr|. (19) 

or 



H 1 + [H o ,S} = 0. (13) 

In Eq.(12), if we discard the higher order terms and only 
keep the second-order term, the effective Hamiltonian can 
be achieved approximately as 

H cS = H Q + ^[H h S]. (14) 

From the Eq.(13) we certainly know how to construct the 
anti-Hermitian operator S, which has the following form 

g =E i ffl)m F WH. (is) 

'—f h, m — £/ n 



H cS = ^2H(n)\n)(n\ : 

n 

H(n) = (hw 3 + ng 2 A_) (a\ + a\) (20) 
Q 2 

In general this is a typical effective Hamiltonian leading 
the two-qubit quantum logic gate. In usual it is obtained 
by adiabatically eliminating the variable data bus with 
various methods |12(. However, in most of previous works, 
this crucial terms 

g 2 A^a(al + a 2 2 ) (21) 
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referred to the ac Stark effect (a dispersive frequency shift 
effects) has been either irrationally ignored or passed over 
in silence. This is unsatisfactory even though we can now 
prove that it can run the logic gate of two qubit system 
in next section. 

IV. A NOVEL DECOHERENCE MECHANISM: 
THE INVERSE STERN-GERLACH EFFECT 

Having gotten the effective Hamiltonian with a perfect 
inter-qubit interaction, we can show how to create the 
quantum entanglement by controllable the coupling be- 
tween photon and qubit. As a data bus, the role of cavity 
field is to introduce extra controllable parameters. In the 
next section we will show the details to implemental an 
ideal two qubit logical gate operations in the decoher- 
ence free subspace (DFS) fS}. However, for those states 
outside the DFS, we can demonstrate a novel decoher- 
ence phenomenon related to the so-called inverse Stern- 
Gerlach|l6| effect from the adiabatic variable separation 
based on the BO approximation [13l |. 

Let us generally consider the adiabatic evolution of two 
identical charge qubits 1 and 2, coupled to a single- mode 
field in the microwave cavity. In the off-resonance case, 
the motion of the qubits does not excite the transitions 
from a cavity mode to another, and then the photon num- 
ber is conserved, i.e., 

[fleff,ata]=0 (22) 

This shows that the total wave function will adiabatically 
keep the factorized structure 

Mt)) = |0„(t)) ® |n) (23) 

during evolution only if the cavity is exactly prepared ini- 
tially in a single number state with definite phonon num- 
ber, namely, a Fock state \n) . In this case the qubit part 
is just governed by the effective Hamiltonian (n\H c $ \n) 
= H(n) and then we can manipulate the qubit system 
according to the n-dependent Hamiltonian to form max- 
imal entanglement. However, if one can not prepare the 
cavity in a single Fock state, the part \<fr(t)) of qubit must 
depend on the different phonon number n and then we 
can not make an exact manipulation for qubit part due 
to this correlation to cavity field. This kind feature of 
quantum adiabatic entanglement is just a novel physical 
source of the quantum decoherence in the process of two 
qubit logical gate operations. 

We remark that this phenomenon is an analog of " in- 
verse Stern-Gerlach effect" in atomic optics, in which 
discrete atomic trajectories are correlated to different 
photon numbers in the cavity. When the atom is non- 
resonant with the cavity modes, there appears a disper- 
sive frequency shift effects affecting both the atomic tran- 
sition and the field mode. It can be interpreted as sin- 
gle atom and single photon index effects. These effects 
lead to various interesting potential applications, which 



have been investigated in cavity QED for atoms, both 
theoretically and experimentally, e.g., the interference 
schemes to measure matter-wave phase shifts produced 
by the non resonant interaction |l6|. This schemes per- 
forms a quantum non-demolition measurement of photon 
numbers in a cavity, at the single photon level. Its ex- 
perimental demonstration is based on the detection of 
Ramsey resonances on circular Rydberg atoms crossing 
a very high Q cavity. For this kind of "inverse Stern 
Gerlach effect", we even presented an extensive gener- 
alization based on the Born-Openheimer approximation 
to analyzed the the adiabatic separation induced quan- 
tum entanglements |l3l |. Thus it defines the adiabatic 
quantum decoherence in general case. We can discuss 
this effect for a solid state based system with two charge 
qubit. 

To see the quantum decoherence due to the general- 
ized "inverse Stern Gerlach effect", we assume the two 
JJ qubits and the cavity field are initially prepared in a 
factorizable state: 

|¥(0)) = |$®|y>) (24) 

where \<fi) is the initial state of the two J J qubits and \cp) 
the state of the field. In general, if the cavity is prepared 
initially in a superposition state of Fock state 

i^} = E c «i n > ( 25 ) 

n 

rather than a single Fock states, the total system will 
evolve according to 

|*(*)> ®l»> (26) 

n 

where 

|0„(t)) = U n (t) |0(O)) . (27) 

The effective evolution matrix U n (t) — exp[— iH {n)t\ is 
governed by ff(n). This "inverse Stern Gerlach effect" re- 
sult from the dependence \<fi n (i)) to different n. 

Let \m) be the single Fock state that we want to pre- 
pare and H(m) be the controlled Hamiltonian. Then we 
can characterize the difference between the real evolution 
and the ideal one p m (t) — \<f> m {t)) ((f> m (t)\, by the fidelity 

F = Tr[ Pm (t) p (<)] = J2 M 2 K<U*)K(*)> I 2 (28) 

n 

where 

p(t) = Tr C (*(t)|) = E l C «| 2 l^n(*)) (*n(*)l 

n 

(29) 

is the reduced density matrix of the two JJ qubits. 

Usually it is difficult to prepare the Fock state |m) 
and we can only use the coherent state \a) with aver- 
age photon number (a|cra |a) = m. Then we assume the 
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junctions are initially in the state |0) 1 |0) 2 , and the ini- 
tial state of the micro- wave field is the coherent state \a). 
Then total system will evolve into 



n 



l*> = e-iW'^|^(t)>| n>: 

n 

= («) 10)! |0> 2 + * Cl («) II)! |1> 2 }, (30) 
where the timc-depcndcnt coefficients arc 



ci (n) = sin (u n t) cos (9 n ) , 
C2 (n) = cos (w„i) — i sin (w„i) sin (6 n ) 




uj n = \\j(2hw i + 2 ff 2 A_n) 2 + (yA+)2, 
sin (0„) =- 



FIG. 2: Fidelity F as a function of time t, with different m, 
which is the mean eigenvalue of the micro-wave field, m = 
0.2(line), m = 0.4(circles) and m = 0.7 (crosses). 



(2^j + 2.g 2 A_n) 2 + (yA+) 2 



Through a simple calculation, we obtain the fidelity 



F = e^l" 1 ! 



{cos 2 (u) n t) cos 2 (u m i) 



+ sin 2 (u n t) sin 2 (0 n ) sin 2 (w TO i) sin 2 (9 m ) 

+ cos 2 (w TO i) sin 2 (uj n t) sin 2 (0 n ) 

+ cos 2 (u) n t) sin 2 (w m t) sin 2 (9 m ) (31) 

+ sin 2 (w„i) cos 2 (6 n ) sin 2 (w TO i) cos 2 (6 m ) 

+ 2 sin (w„f) cos (u) n t) cos (#„) sin (u> m t) cos (0 m ) cos (w m i) 

+ 2 sin 2 (w„i) sin 2 (u m t) cos (0„) cos (0 m ) sin (9 m ) sin (6>„)}.40(circles) and t = 70 (crosses) 
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FIG. 3: Fidelity F as a function of the mean eigenvalue of the 
micro- wave field m, with different fixed time t = 13(line),t = 



FIG. 2 and FIG. 3 illustrate that the fidelity decays 
sharply with the value of to, namely, coherent state |a) 
leads to big deviation of p a and p m with big a. This is 
because that coherent state \a) is just in Fock state \m) 
with the probability P m , 



The real evolution process governed by the effective 
adiabatic Hamiltonian H(k) is 



(*)> 



|00) 



4 (fc)|oo) + ici (fc)|ii», 



p m\ a = m = \(m \a) 



= e"H 



2 TO 



to! 



When to — ► oo, P TO — > quickly. 

The above discussion shows us that, when we prepare 
the controllable cavity field in different initial states, one 
can get different entangled states for the two qubit. This 
motivates us to explore the possibility to realize the per- 
fect logic gate operation by initially preparing the micro- 
wave field in coherent state and Fock state. Let us con- 
sider the above mentioned problem in the following. 

We aim to get a standard Bell state \(j> + ) 

|0+) = i=(|OO) + |ll» 

from the state |0) 1 |0) 2 , by preparing the cavity field ini- 
tially in Fock state |ft). The evolution from 1 0) x 1 0} 2 to 
\4> + ) naturally realize a perfect ideal logic gate. 



where the Hamiltonian H(k) corresponds to the Fock 
state | ft) for fixed ft. We can use the square of the norm 
of the inner product 



/ fe H<cw>r 

to characterize the difference between the ideal state 
and the real state \<f> k ) 

|2 



: 0(1 



(t)\4> + )\ 



(32) 



cos 2 (u) m t) + sin 2 (u m t) [cos (9 rn ) + sin (6 m )Y 



The above equation shows that f m is a periodic function 
of time t. 

FIG. 4 shows that f m decays with the average photon 
number to and the maximum value of f m can not ap- 
proach 1. Based on this result, we can not construct an 
ideal logic gate in this system. 
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FIG. 4: The vertical axis represent the function f m , the hori- 
zonal axis represent time t, m is the eigenvalue of the Fock 
state, m = (line), m = lO(crosses) and m = 20 (circles) 




FIG. 5: The vertical axis represent the function f a , the hori- 
zonal axis represent time t, and a is the eigenvalue of the 
coherent state, a = 0.1 (circles), a = 1.1 (crosses) and a — 5 
(line) 



Now we consider another case that the micro- wave field 
is prepared in coherent state \a) initially, and the junc- 
tions is prepared initially in the state |0) 1 |0) 2 . By a 
simple calculation, we get the final state depicted by the 
reduced density matrix(RDM) 



Pa 



tr mw {\V) (*|} 

l |2fe 



It is shown from the FIG. 4 and FIG. 5, whatever the 
state of the microwave field is prepared in, the ideal logic 
gate operation can not be realized in this system. But 
this does not means that we can not obtain any maximal 
entangled state in this way. We will discuss this problem 
in next section. 



V. CREATING MAXIMAL ENTANGLEMENT 
IN THE DECOHERENCE FREE SUBSPACE 

From the discussions in the above section, we find that 
the maximally entangled state can not be obtained only 

1^ |1) 2 , even though one can 



= e-l«l 



k\k\ 



(33) 



from the state |0) 1 |0) 2 and 
prepared external controlled microwave cavity field in an 
arbitary state. While the other two states |0) 1 |1) 2 and 
|0) 2 can span a decoherence-free subspace(DFS)iHT 1 , 
it means that any superposition of the state |0) 1 |1) 2 and 
|1) 1 |0) 2 can evolve into this kind of DFS. Easily seen in 
Eq.(19), the effective interaction between the cavity and 
qubits g 2 A_a^a (<r\ + cr 2 ) vanishes in the DFS and can 
not distinguish between any two states in this DFS. So 
we conclude that there is not a "which- way detection" 
to determine the "paths" in this case, i.e., there is not 
decoherence appearing in DFS. 

Let us use a special example to demostrate the above 
observation. When junctions are prepared initially in the 
state |0) 1 |1) 2 and the manipulative field prepared in the 
Fock state |n), then the evolution of the total system will 
evolve into 



\M\n) = e~ im 10)^1)^) 

= cos^A.^lO^li^ln) 

-isin( 5 2 A_t) |1> 1 |0) 2 \n) (35) 

It is obvious that, when t = —^^2^y^ ; thc two qubit 
system reaches a maximally entangled state 

|^oi> = ^(|0> 1 |l> 2 -i|l> 1 |0> 2 ). (36) 

By the same way, we can obtain the other maximal en- 
tangled state 



We explicitly calculate the function 

f a = \(<t> + \p a \t + )\ 



= l e -M s 



i2fc 



km 



| cos 2 {(jj k t) 



+ sin 2 (ukt) [cos(6» fe ) +sin(6» fc )]' 



(34) 



FIG. 5 displays the evolution of the function f a calcu- 
lated from eq. (34) for the qubits with coherent microwave 
fields. When the eigenvalue a — > oo, the fidelity F — » . 



1 

V2 



IV'io) — (10)! |1) 2 +*!!>! |0) 2 ) 



(37) 



when we take the t = ^^^^'^ ■ Both IV'oi) an( i IV'io) 
are independent of the controllable micro- wave field, they 
belong to the DFS 2U 1 . The basis 10)! |0) 2 and |1) 1 |1) 2 
span the other subspace of the Hilbert space of the JJ 
qubits(2U), we denote W]_. Thus we have 2U = W 1 ®^. 
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VI. OVERALL QUALITY OF CREATED 
ENTANGLEMENTS 

In above section we have discussed that when junc- 
tions are prepared initially in the state |0) 1 |0) 2 or state 
|1) 2 , we can not obtain maximally entangled state of 
the junctions with any controllable microwave field. But 
we can study the entanglement of these states evolved 
from the initial state jO), |0) 2 or state |1) 2 . 

As well-known, for a bipartite system, composing of 
two subsystems A and B, the bipartite entanglement can 
be measured by its concurrence 1 1 7j which is defined by 

C(p) = max(0, Ai - A 2 - A 3 - A 4 ) (38) 

where Ai, A 2 , A3, A4 is the square root of non-Hcrmitian 
matrix R in decreasing order, and 

R = p(a v ®a v )p*(a y 1 ®a y 1 ). (39) 

We consider that when the initial state of junctions is 
|0) 1 1 0)2 and the controllable microwave field is prepared 
in coherent state. Then, through a simple calculation, 
we obtain the concurrence of the states of JJ qubits sub- 
system Eq.(33) is 




FIG. 6: The vertical axis represent concurrence C, the hori- 
zonal axis represent time t, a — 0.1 (light ashen line), a = 1.1 
(ashen line) and a = 3 (black line). 

By a simple calculation, we get the RDM of the J J qubits 
Pn = ^£^ nPE ^n){^ n \ ■ (43) 

n 

By the same way, we calculate the concurrence of this 
state in the following 




A = e -l«l 2 ^^L_| cos(c , fe t)_ iS i n ( Wfci ) sin ( 0fc )|2 

k 

B = er^Y^^n 2 (uj k t) cos 2 (9 k ) 

k 
k 

■ {cos (&k) cos (u>kt) — i sin (tUkt) sin (9k)} • (40) 



In FIG. 6, the concurrences of the qubits are plotted for 
different values of a. It is seen that, with the increasing 
the eigenvalue of the controllable microwave field, the 
concurrence decrease sharply. 

We study another case, when the initial state of the 
controllable micro-wave field is thermal state 

Pmw = ^J2 e " f3E l n ) ( n \ ^E e ~"^ ( 41 ) 

n n 

the state of the total system evolute into 

P=^Y, e ~ n0E ^n)\n){n\{i> n \ ■ (42) 

n 




k 



B = -Y,e~ k0E sin 2 (^kt) cos 2 (0 k ) 

k 

D = I^Te-^sinKOsin^) 

• {cos (Ok) cos (uiki) — i sin (uikt) sin (6k)} 
In FIG. 7, the concurrence of the state pjj is periodic 

function of time t and the concurrence C — > the maximal 
value, when (3E — > 00. This is because f3E — * 00, the 
thermal state p mw — > |0) (0|. 

FIG. 6 and FIG. 7 illustrate that any types of the con- 
trollable micro-wave field can not increase the entangle- 
ment of the JJ qubits when its initial state is superposi- 
tion of 1 00) and |11), and the maximal entanglement is 
much smaller than 1. Only the initial state of the control- 
lable micro-wave field is vacuum state, the entanglement 
can reach the maximal value. 



VII. CONCLUSION 

In summary, we study the protocols which can create 
maximally entangled states between two qubit coupled 
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FIG. 7: The vertical axis represent concurrence C, the hori- 
zonal axis represent time t, with different parameter/3 E of the 
thermal states. (5E = 0.7 (light ashen line), (5E = 2 (ashen 
line) and /3E = 6 (black line) 



to a controllable microwave field in a cavity. In order to 



obtain the analytic study for this decoherence problem, 
we generalized Frohlich transformation to re-derive the 
effective Hamiltonian of these system, which is equivalent 
to that obtained from the adriatic elimination approach. 
Because of nontrivial decoherence, we can not construct 
an ideal logic gate by this system. But we can construct 
a decoherence-free subspace of two-dimension to against 
this adiabatic decoherence in this system. 
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